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Abstract 

Travelling waves of densities of binary fluid mixtures are investigated near a critical 
point. The free energy is considered in a non-local form taking account of the density 
gradients. The equations of motions are applied to a universal form of the free 
energy near critical conditions and can be integrated by a rescaling process where 
the binary mixture is similar to a single fluid. Nevertheless, density solution profiles 
obtained are not necessarily monotonic. As indicated in Appendix, the results might 
be extended to other topics like finance or biology. 

Key words: Binary fluid mixtures; critical points; travelling waves; rescaling 
process. 
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1 Introduction 



In physical chemistry, thermodynamics and condensed matter physics, a criti- 
cal point specifies the conditions (temperature, pressure and concentration) at 
which distinct phases do not exist [6,9,29]. There are multiple types of critical 
points such as vapor-liquid or liquid-liquid critical point. A single fluid has a 
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unique critical point associated with given temperature, pressure and density. 
For binary mixtures of fluids, in the space of temperature, pressure, concentra- 
tion, critical points are represented by a curve in a convenient domain [22]; to 
each temperature we can associate a critical pressure and two critical densities 
corresponding to the mixture components [T3|36] . 




An important thermodynamical potential is related to the mixture volume 
free energy [32 33 34J. At a given temperature, the volume free energy is as- 
sociated with the spinodal curve connecting the two different phases of the 
binary mixture. Due to conditions of equilibrium of phases, it is possible to 
form a general expansion of the free energy near a critical point. This form 
is known in the literature by means of physical chemistry considerations [37] 
and is the form we use in our calculations. 

By calculations in molecular theories, the densities of the components fluc- 
tuate near a critical point [24]. In the following, we use a continuous model 
to investigate how the average variations of densities are related to molecular 
interactions. Two assumptions are explicit |4]I18|38] : 

i) The component densities are assumed to be smooth functions of the dis- 
tance from an interface layer which is assumed to be flat on the scale of 
molecular sizes. The correlation lengths are assumed to be greater than inter- 
molecular distances [2T|23] ; this is the case when at a given temperature T 
the parameters are close to the ones of a critical state [35] . 

ii) The binary mixture is considered in the framework of a mean-field theory. 
This means, in particular, that the free energy of the mixture is a classical 
so-called "gradient square functional". This kind of Landau- Ginzburg model 
consisting of a quadratic form of the density gradients comes from Maxwell 
and van der Walls original ideas [27|31|42|T44] . At given critical conditions, the 
coefficients of the quadratic form are constant. 

This point of view that, in non-homogeneous regions, the mixture may be 
treated as bulk phase with a local free-energy density and an additional con- 
tribution arising from the non-uniformity which may be approximated by a 
gradient expansion truncated at the second order is most likely to be suc- 
cessful and perhaps even quantitatively accurate near a critical point [37] . 
The approximation of mean field theory does provide a good understanding 
and allows one to explicitly calculate the magnitude of the coefficients of the 
model. These non-linear equations are able to represent interface layer and 
bulks and consequently allow to build a complete theory of the mixtures in 
non-homogeneous domains in dynamics. 

In Section 2 we recall the equations of motion in a pure mechanical process 
obtained through the Hamilton variational principle. 

Section 3 is devoted to travelling waves without dissipation. Due to the fact 
that the equations are Galilean invariant, the case of equilibrium and the case 
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of motion are analyzed together. 



In Section 4 by means of a rescaling process taking the vicinity of a critical 
point into account, we integrate the equation for equilibrium as well as for 
motions with dissipation. 

Two appendices present the motion equations and the mathematical reason 
of the choice of the free energy form near a critical point of a binary mixture 
of fluids obtained by a new method issued from differential geometry. 



2 Isothermal motion of a binary fluid mixture near a critical point 

We study a mixture of two fluids by a mechanical process. No assumption 
has to be done about composition or miscibility. The motion of a two-fluid 
continuum can be represented with two surjective differentiable mappings (see 



where subscripts 1 and 2 are associated with each constituent of the mixture. 
Term z = (t, x) denotes Euler variables in space-time W and terms Xi and 
X2 denote the Lagrange variables of constituents in reference spaces T> \ and 
V o2 respectively. 



Fig. 1) 




z -)■ Xi = $i(z) and z -» X 2 = $2(2) 



D, 



'a I 
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Fig. 1. General representation of a two-fluid continuum motion. 
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In the pure mechanical case, the Lagrangian density of the mixture is 

where Vi and v 2 denote the velocity vectors of each constituent, pi and p 2 are 
the densities, Oi and f2 2 are the external force potentials depending only on 
z = (t, x) and e is the volume energy [2~|TT4"] . 

The expression of the Lagrangian is in a general form. In fact dissipative 
phenomena imply that vi is almost equal to v 2 ; it is the reason why we do 
not take account of some kinetic energy associated with the relative velocity 
of the components which is of smaller order (at least of order 2) and will 
be negligible in travelling wave behavior [5125128] . Because of the interaction 
between the constituents, the volume energy e is not the sum of the energies of 
each constituent of the mixture, like for Euler mixtures of fluids. The mixture 
is assumed not to be chemically reacting. Conservations of masses require 

Pidet F< = poi (Xi), (1) 

where subscript i belongs to {1,2}. At t fixed, the deformation gradient <9xj/<9Xj 
associated with is denoted by Fj and p oi is the reference specific mass in 

Eqs ([T]) are equivalent to the Eulerian form 

^ + divp jVj = 0. (2) 
at 

The volume energy e is given by the behavior of the mixture [71111112] . In our 
mechanical case, for an energy depending on gradients of densities, the volume 
energy is 

e = e(pi, p 2 , gradpi, gradp 2 ). 

The potential 

de d ( de \ 
^ % dpi <9x 7 \dpi n ) 
defines the specific free enthalpy or chemical potential of the constituent % of 
the mixture [16j. Subscript 7 corresponds to the spatial derivatives associated 
with gradient terms. Usually, summation is made on repeated subscript 7. In 
practice, we consider a quadratic form with constant coefficients Ci,C 2 ,-D 

Q = Ci (gradpi) 2 + 2£>gradpi gradp 2 + C 2 (gradp 2 ) 2 

such that 

e = 9o(pi,p2) + 7; Qi ( 3 ) 
where g (pi, P2) is the value of the volume energy of the homogeneous bulks. 

To obtain the equations of motions, we use a variational principle whose orig- 
inal feature is to choice variations in reference spaces (Fig. 1). They are as- 



4 



sociated with a two-parameter family of virtual motions of the mixture (see 
Appendix A). 

The equation of the motion of each constituent of the mixture writes ([12] and 
the references therein) 

^ + grad (jn + SU) = 0, % = {1,2}, (4) 

where aj denotes the acceleration of the component i (2=1,2). In applications, 
the motions are supposed to be isothermal (T denotes the common temper- 
ature value of the two components) and correspond to strong heat exchange 
between components. In thermodynamics, this case corresponds to a function 
9o{P\i P2) as the volume free energy of the homogeneous mixture at tempera- 
ture T. 

In our model, the equations of motion (TjJ yield 

ai = grad {d Ap 1 + D Ap 2 - g , n } 

(5) 

a 2 = grad {D Ap 1 + C 2 Ap 2 - 9 o, P2 } ■ 

Taking Eq. (J2]) into account, we can note that the two equations of system (jlj) 
are equivalent to the system 

d v 

-^—^ + div (pM ® Vj) + grad {pi = & gradp i; i = {1, 2} 
at 

for which equations of components are not in divergence form, but the summa- 
tion of the two equations and the fact that Y%=i A*i grad pi = grad e allow to 
obtain the total motion of the mixture in divergence form [MfTU] . Therefore, 
while the global momentum equation represents a balance law, individually, 
the equations of system §2§ do not. 



3 Travelling waves for an isothermal fluid mixture in one-dimensional 

case 

3. 1 The system of equations for travelling waves 

An interface in a two-phase mixture is generally schematized by a surface 
without thickness. Far from critical conditions, this layer is of nanometre size 
and density and entropy gradients are very large. This is not the case in the 
immediate vicinity of critical points where fluctuations of densities strongly 
diverge. Our goal is to schematize the average of these fluctuations by means 
of travelling waves in a continuous model using an energy in form (J3]). 
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To obtain one- dimensional travelling waves in isothermal case when we neglect 
diffusion and viscosity, we are looking for solutions only depending on the 

f = x - st, 

where s denotes the wave velocity. Equations (J2J) of conservation of masses 
write 

dpi d(piVi) = Q 
dt dx 

where Vi denotes the one-dimensional velocity of component i. Equivalently 
-7?{Pi{Vi ~ s)} = ^=^> Vi = — + s, (6) 

GEcJ Pi 

where d{, i = {1,2}, are constant along the motion and have the dimension 
of a mass flow. The acceleration of mixture component i is 



dvi dvi 
ox dt 

and system (jSJ) yields 



d_ ( (vj - sf 
dU 2 



ld_ (of 

2dx [pi 



Id 2 

C lP '[ + Dpl = pi + -^ + h 
z p 1 



Id 2 

Dp'[ + C2P2 =P2 + ~^ + k 2 
* P2 



(7) 



where k\ and k 2 are two constants, superscript T denotes the transposition 
and 



X 



Pi \ D Ci D 1 d\ 1 d\ , 

, R = , Y(pi, p 2 ) = g {pi, Pi) -7; o —+hPi+k 2 P2- 

P 2 \D C 2 2 Pl 2 p 2 



Here ' = d/dC, denotes the derivation along the ray. Let us note that the sum 
of equations of system (jTJ), respectively multiplied by p[ and p' 2 , yields the first 
integral 



Id 2 \d\ (\ 



9o(pup2) + h pi + k 2 p 2 - J -± - J ^ - ( J d p'( + D Pl p 2 + J C 2 p'( ) = fc 3 , 

1 p\ 1 p 2 \1 2 / 

where fc 3 is a constant. 

Consequently, in one-dimensional system, the profiles of mixture densities are 
associated with a mechanical system having an equivalent kinetic energy 

1 „ In _ / / 1 _ In 

~C X Pl 2 + D Pl p 2 + -C 2 p 2 
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and an equivalent potential of forces 

Id 2 Id 2 

/ \ 7 7 1 2 

(Pi, P2) - h pi - k 2 p 2 + - — + - — . 

Z pi I p 2 

It is possible to use all the tools of analytical dynamics as Lagrange or Hamil- 
ton equations, symplectic geometry or Maupertuis methods. 



In each bulk, density gradients are null. Eliminating constants ki,i = {1,2}, 
dynamical conditions through the interfacial layer yield 



[g, P i(Pi, Pi)} = 

[g, P2 (pi, P2)} = 
[g- Pig, Pl - p2g, P2 ] 



where [u] = Up — u a denotes the difference of values of u between the two 
mixture bulks a and /3 and 



9 = go 



1 eg 
27i 



14 

2 p 2 



In case of equilibrium at a given temperature T, (di and d 2 are null), the 
minimum of the total free energy, with a given total mass for each constituent, 
yields conditions 



Mpi,p 2 )] 

[P{PUP2)\- 




= 
. 



(9) 



where 

P(Pl,P2) = Pc + PlPliPl, P2) + P2P2(PI, P2) -go(pl,p2) 

denotes the total thermodynamical pressure of the binary mixture. The con- 
stant P c , corresponding to the fact that g D is defined to an additive constant, 
can be defined as the critical pressure term for a value of g Q null at the critical 
point. 

Conditions of system (jSJ) (or system (jUj)) express that the bulks of Gibbs sur- 
face (£) associated with the dynamical volume free energy g (or static volume 
free energy g Q ) correspond to the contact points with a bitangent plane (Fig. 

1 d 2 1 d 2 



2). By adding the term — - 



— to g , the study of non-dissipative trav- 
2 pi 2 p 2 

elling waves of density components of a mixture turns back to an equivalent 
equilibrium problem. 
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3 (So) 




i 



Fig. 2. In case of bulk equilibrium, the third coordinate of the point of intersec- 
tion of the bi-tangent plane with the (/-axis is —P{p\, p^j. In dynamical case, the 
intersection corresponds to the opposite of a dynamical pressure (and defined to 
an additive constant P c ). Contact points M a and Mp of the bitangent plane with 
surface (£) generate a curve (T). Along this spinodal curve (T), the straight line 
(-D) connecting the contact points M a and Mr associated with the two mixture 
bulks is the characteristic line of the plane which generates a developable surface. 
Critical point C corresponds to the collapse of the two contact points at the given 
temperature T. 

4 Integration of the equations of travelling waves 

In Appendix B, we propose an expansion form of the volume energy g for the 
binary mixture, near its critical point at temperature T, in the form 



where A and B are physical constants, q\ and qi are convenient linear combi- 
nations of pi — p\ and p 2 — P2 with pi and p c 2 denoting the density values of 
the two constituents at the critical point (for the sake of simplicity, we denote 
the volume free energy as g Q in all systems of variables). With the change of 
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variables, 



pi - pi 

p2 - P% 



(10) 



and 



Y(r u r 2 ) = g (n,r 2 ) - - 



1 d\ 1 d\ 



2 ri + pi 2 r 2 + p 2 



- + h ( n + pi) + k 2 (r 2 + pi) 



system (0) can be written 



Rr"= -7T- 



<9r 



We denote 



where scalars a, b, c, c? are depending on the thermo-mechanical properties of 
the mixture near the point C . By the new change of variables, 



r = P q , with q 



'a 



(12) 



and 



^(9i,92) =^0(91,92) 



1 



1 



d\ 



2 agi + 6g 2 + Pi 2 cgi + G?g 2 + p 2 
+ki {aqi + bq 2 + p\) + £; 2 {cqi + dq 2 + p%) 



the system (fTTT) writes in the form 



Rq" 



'fr) 



with 



R 




or 



CVg? + D q'i = AAB 2 (B 2 qf - qi q 2 ) + k x 
D q'l + C 2 q' 2 ' = 2A (-B 2 qf + 2g 2 ) + k 2 , 



(13) 



where ki = aki + ck 2 and k 2 = bk\ + dk 2 . Due to the fact R is symmetric 
positive definite and P is not singular, the same holds for R = P T RP. 
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4-1 Case of equilibrium 



4-1.1 The rescaling process 

The values of constants k\ and k 2 must correspond to the singular points of 
the differential system (J7|) when d{ = 0, (i = 1,2). These singular points are 
associated with the two bulks of the binary mixture. Consequently, k\ = k% = 
and k 2 = k 2 = —t 2 (see Appendix B, for the meaning of "distance" to 

the critical point). 

We consider the following hypothesis (Hi) and (i?2) : 

• (Hi) We assume that gi and q 2 slowly vary as a function of the one- 
dimensional parameter x. The hypothesis can be expressed by a change of 
variable z = e x, where e is an adimensional small parameter (e 1). With 
the following change of variables: 

qi(x) = Qi(z) and q 2 (x) = Q 2 (z) , 

system ( fl3|) yields 

e 2 (C1Q1 + DQ 2 )= AAB 2 (B 2 Q\ - QiQ 2 ) , 
e 2 (p Qi + C 2 Q 2 ) = 2A (~B 2 Q 2 + 2Q 2 ) - r 2 , 

where the dot denotes the derivation with respect to z. The solution must be a 
phase transition connecting singular points M a (r/B,T 2 ) and M a (—r/B,T 2 ). 

• (H 2 ) The density profiles are of small amplitude with respect to the critical 
densities. Consequently, we look for the solutions in the form 

Qi(z) = e^yi(z) and Q 2 (z) = e n *y 2 (z) , 

where ni and n 2 are to real positive constants. System (fl4"l) yields 

e 2 (Cie ni y\ + D e n2 y 2 ) = AAB 2 (B 2 e 3ni yf - e^ +n2 yiy 2 ) , / 

;~ ~ ( ( 15 ) 

e 2 (D e ni y\ + C 2 e n *y 2 ) = 2 A (-B 2 e 2ni y 2 + 2e n '*y 2 ) - r 2 . 

The form of the volume energy near the critical point must be conserved by 
affinity This assumption corresponds to the universality of the form of the 
energy near a critical point |10| . Consequently in this rescaling process, 

2m = n 2 and r = e ni E , 

where E is a positive constant with physical dimension M 1 / 2 ^ 3 / 2 : . System 



10 



( TT5|) yields the equality of same order terms, 

' s 2 +^Cm = e 3n HAB 2 (B 2 yf - ym ) , 

~ (16) 
£ 2+m Dyi = e 2ni 2A {-B 2 y\ + 2 y 2 - E 2 ) . 

If we assume 2 + ni < 2n\, Eq. f fl6|) 9 reduces to ?/i = 0, which is not phys- 
ically possible; If we assume 2 + rii = 2ni, Eq. f JT6|) 1 reduces to j/i = 0. 
Then, to obtain non-trivial solutions, we must assume that 2 + rti > 2 ri\ and 
2A (-B 2 y 2 + 2y 2 - E 2 ) = 0; consequently, 

Ifc = ~ V + E 2 ) . 

Now, if we compare e 2+ni and e 3 ™ 1 , we must assume that n\ = 1. System (Fl6l) 
leads to the reseated system 

' C lV \ = 2AB 2 (B 2 yf-E 2 yi ) 
\ V2 = \(B 2 y 2 l+ E 2 ) 

which is equivalent to 

jt^^-sv+fa) (17) 

where is a convenient constant of integration. 



4-1.2 Integration of system (17) 

The first equation of system (1171) has several possible solutions. In this case we 
look for phase transition, the solution must be monotonic and must go from an 
equilibrium point to the other one. The simplest possibility is = E 4 /(2B 2 ). 
The integration of Eq. (TTTj) i yields 



yiK z ) 



E 
B 





2AB 2 \ 


— - — X 







Consequently, with the change of orientation of the x-axis, we get 



0i W 



Q2{X) 





2AB 2 \ 




— - — X 




c x ) 



T 

B 

— { I 4- tanh" I r, 
2 1 \ \ Ci 



2AB 2 
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Due to the changes of variables ({TO]) . (fl2|) . we obtain the general form of the 
phase transitions 



Pi(x) = pi + a\ tanh ( — ) + b\ \ 1 + tanh 



x 



p 2 (x) = p c 2 + a 2 tanh (^) + ^ j 1 + tanh 2 (jr\\ ■ 
Depending on the two cases of boundary conditions, we obtain 
I s * Case: When lim pAx) = p" and lim pAx) = pf with i = 1,2 , then 

x— oo x— >+oo 

ai = l(pP-pf), b l = \{pt-2pl + p^. 
2 nd Case: When lim pAx) = p^ and lim pAx) = pi with i — 1,2 , then 



Let us note that the constants and frj depend on the thermodynamical 
behavior of the mixture and they may be positive or negative. 



4-1.3 The different cases 

The variation of pi(x) and p2{x) in system ffl8l) yields the different cases. They 
are summarized on Fig. 3 as domains delimited by the bisectrices of axes 
and hi and we represent the different possible variations of each component 
on Fig. 4. We it is known by other methods in physical chemistry that 

the density profiles in interfacial densities are not necessary homogeneous for 
mixtures of fluids [37f38] . 



4-2 Case of travelling waves 



In the case of motion, we check with travelling waves submitted to viscosity; 
in this case, system (Q can be modified to take account of the viscosity in 
each component of the mixture 

( 1 

ai = grad {Ci Ap 1 + D Ap 2 - g } H div <j x 

Pi 

< 

a 2 = grad {D Ap 1 + C 2 Ap 2 - g } + — div a 2 , 
k - p 2 
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Fig. 3. The different possible cases for each mixture component are presented on 
the graph associated with the aj-axis and 6j-axis. Points on axis a,, where pf = p^ 
are forbidden; this will not be the case in dynamics. 







Fig. 4. Sketch of the different possible forms of phase transition for each mixture 
component are presented on the graph with densities pi versus £. 



where Oi{i = 1,2} is the viscous stress tensor associated with component i. 
Consequently we obtain in variables q\ and q 2 , the system of equations of mo- 
tions in one-dimensional system with cr, = v", where z/j > is the coefficient 
of viscosity of the component i [39]. The system can be written by using one 
dimensional derivation in diver,, by integration as in section 4.1 and with the 
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fact that pi 



Pi, 



C iq '{ + D q% = AAB 2 (b' 



lac?? 



2 (p9) 2 



1 c d 2 



(PIT 

D q'l + C 2 q' 2 ' = 2A (-B 2 q\ 



qiq2 

I Cd; 



(ft. 



— ?2 + «1 
'lJ 



l_6_dj_ 

+ 2(^ 



1 d d 2 

+ - 



2 (Pi) 2 (p!) 



'9i + 

bdyvy , dd 2 v 2 f j 



(19) 



4-2.1 The rescaling process 

We can make an analogous rescaling process as in Section 4.1.1. The differ- 
ence only comes from scalars d\ and d 2 - They must be of small amplitude; 
consequently, we look for constants in the form 

D 1 = e n3 di and D 2 = e" 4 d 2 , 

where constants and are positive. To obtain convenient solutions we 
must have 



di 



3ni 

Die 2 , U 



1,2), k x = e 3ni K u k 2 = e 2ni K 2 , n x = 1 



Then, the terms — — \q\ = — ^—^ e( 3+2 ^yi and 2 ! ol = — 2 -^ 2 - e( 4+2 )y 2 

(pi) 3 1 (Pi) 3 (Plf 2 (Pi) 3 

which are of order greater than the other terms in system ( fl9l) . are negligible. 
Finally, system ( f!9l) yields 



Ciy? = AS 2 (fi V + ^y 2 + K m + 



(20) 



where K 3 = 
integration. 



a D 2 



+ 



c 



AB* \ 2 (pi) 2 2 (p 



+ Ki and is a new constant of 



Integration of system 

A qualitative study of Eq. ( I20l) i can be considered. Constant A is strictly 
positive. The discussion comes from 

y\ = f(yi) -g(yi), 
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where 



f(yi 



ALV ( 2 t K 2 V 



Ci 



(21) 



represents a quartic and 



AA 2 B 2 



(22) 



represents a straight line. If K 2 > 0, f(vi) is a convex quartic and only periodic 
solutions are possible. Consequently K 2 < 0. 

Two interesting cases can be considered: 

a) The straight line (|22|) is bitangent to the quartic (12~TT) . The solution of Eq. 
(I20l) i is a travelling wave of phase transition in the same form than the two 
first curves of Fig. 4. We must have 



K 3 = 0, K 4 



AA 2 B 2 



and consequently K\ 



aD\ 



cD 2 



2 (pi) z(pi 



\2 • 



Then, 



9i(0 



x 
B 



tanh 



2AB 2 



,/of-i = _^l + tanh 2 |x^-~- 



with x = e\ 



-K 2 
2A 



We note that x plays the role of r in the static case. The solutions in variables 
pi and p 2 are in the same form than in Section 4.1.2 and the same sort of 
density profiles can be observed 



pi(0 = p! 



^ tanh U 



,2, £ 



6x 2 r 

— U + tanh^ 
2 1 V L x 



r cx , / £ \ dx 2 



1 + tanh 2 (-^- 



with L x 



X 



2A5 2 



We have seen that to obtain travelling waves near a critical point, the viscosity 
terms have to be negligible with respect to the others. This result is in accor- 
dance with dispersive waves as the ones obtained by using the Korteweg-de 
Vries equation [43J. 
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5 Conclusion 



The equations of binary mixtures of fluids are considered in the one-dimensional 
case. By using an asymptotic expansion of the volume free energy near a crit- 
ical point and by taking the density gradients into account, we get, in contin- 
uum mechanics, a dynamical system representing the equations of travelling 
waves for the two-fluid components. A rescaling process associated with the 
universality of the energy form near a critical point allows to obtain different 
sorts of travelling waves of density as well as in conservative and dissipative 
cases. We can discuss different possible phase transitions according to the val- 
ues of phases densities of components and the characteristic of the mixture. 

Now, we consider some physical comments. For given densities, at time t — 0, 
the velocities of the mixture components are deduced from the mass balance 
equations ([6]). Then, an arbitrary constant s is introduced in the velocity fields: 



The initial conditions yield this arbitrary velocity s. 

First case: interface propagation - By crossing the interface (interface means 
the domain where the densities drastically change), the volume changes. This 
phenomenon cannot occur in a closed tube with constant volume, but only in 
a tube with one closed end. It is similar to a phenomenon of vaporization or 
condensation. 

One imagines that the other side of the tube is closed with a piston whose 
displacement may be imposed. The component velocities at the fixed end are 
zero in bulk /?; this determines the constant s thanks to the condition 



Due to the fact that p[ ~ p\ and — P21 we S e ^ 

s= d -l = d -l (23) 

Pi P2 

and the values of di and d 2 are not independent. 

Second case: solitary waves - The volume of the interface only moves in one 
bulk phase and remains constant. Such a wave may move in a closed tube as 




and at any time t the relation leads to 




di G?2 
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Natterer's one [6]. At the ends which are assumed far from the interface, the 
velocities in the bulk phase are zero and we still obtain Eq. (|23|) such that the 
values of d\ and d 2 are not independent. 



6 Appendix A: Proof of system ([5]) 



Hamilton's principle - variational form of the principle of virtual powers - 
allows us to derive the equations of motions. The variations of motions of 
particles are deduced from the functions 

X^^zjA), IE {1,2} 

for which are two parameters defined in a neighborhood of zero; they are 
associated with two families of virtual motions of the mixture. The real mo- 
tions correspond to = : ^(z ; 0) = $j(z). 

Virtual material displacements associated with any variation of real motions 
can be written as [TB] 



- — - 

Such a variation is dual with Serrin's |40j. This has been studied in [16] and 
corresponds to the natural variation of the motion in a Lagrangian represen- 
tation. 

Let us denote & = 5;Xj ; we deduce the two relations in Lagrangian variables 
[81I30B0] 

r 9 & r i. c . Pi d Poi , 

ot poi <9Xj 

where <9/<9Xj is the linear form associated with the gradient and div D j is the 
divergence operator on T> oi . Assuming that terms on the edge of W are zero, 
we obtain 



d ( de 



5iPi > dvdt 



or 



$ia= / p i < i?i div oi & H --^^ - vf Fj— - 1 > cfo oi dt 



where i?j = ^v^ — pj — f^. But, 



div oi (p oi i?i^) = p oi i?idiv oi & + -RiT^p 6 + P oi ~^r &■ 
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So, the terms given by integration on the edge of T> Q i being zero, 

„t 2 p ( QR Q ~) 



Finally, for each constituent, 

a (v < F '> = Ml 

and by using the relation 
we obtain 

7, 1 dvf dRi rp d . „ . 

This leads to the vectorial form 

+ grad(/ij + Oj) = 0, i = {1, 2}. 



7 Appendix B: The volume free energy of a binary mixture near 
critical conditions 



In this Appendix, thanks to a direct method of differential geometry, we ob- 
tain the form of the free energy of a binary mixture near a critical point. In 
physical chemistry, the form is presented thanks to the knowledge of the ther- 
modynamic coordinates of the bulks when we are near a critical point |37j . 
In this new presentation, we take only account of conditions of equilibrium 
of bulks when the critical point marks the limit of their coexistence. Such a 
method implies general physical consequences as in one component fluid when 
the two variables are densities of matter and entropy but might also be ex- 
tended in other topics when different equilibrium states collapse only in one 
state as in finance or biology. 

The volume free energy of binary mixture at temperature T is taken in the 
form: 

where p\,p% are the densities of the two components. The following relations 
between the chemical potentials ^1,^2 of the two constituents and g hold 

01 = (9o)' , p>i = (g )' 
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and the thermodynamical pressure is given by 

P = ~9o + PlUl + P2/i2 + Pc = ~9o + Pi (9o)' pi + P2 (9o)' p2 + Pc- 

As the function g is twice differentiable, the geometrical representation of g Q 
by a surface X admits a tangent plane at any point. The normal vector to the 
plane has the direction coefficients (-/ii -/i 2 , 1). 

The phase equilibrium of two bulks a, (3 is described by the conditions 

T(a)=T((3) , ^(a) = , to(a) = M , P(a) = P(0). 

It is noticeable that, in this case, the third coordinate of the intersection point 
of the bitangent plane with the third axis is —P(pi,p 2 ,T) + P c (see Fig. 2). 
The contact points M a , Mp of the bitangent plane to £ generate two curves 
(r a ), (Tp) respectively. The curve (r) = (r a ) U (Fp) is the spinodal curve of 
the phase transition. The case with only one contact point corresponds to a 
mixture with only one bulk while two contact points with a bitangent plane 
correspond to a two-bulk equilibrium mixture. 

When, at a given temperature T, M a tends to Mp, the limit point C is a 
critical point corresponding to the fact that two bulks collapse into one bulk. 
The study of the mixture equilibrium near the critical point C, at a given 
temperature T, is closely related to the shape of surface S in the vicinity of 
C. So it is interesting to prove some theorems giving us information about the 
geometry of the surface £ in the vicinity of C. 

Theorem. The tangent plane to the surface £ at point C is the osculatory 
plane of (T). 

Proof 

Let OAT = F(t) the equation of curve (r) so that F(oj= ot. The tangent 
plane at M a or Mp is 

M a Mp = Fj^5 - F(o5 - F0 + F(o5 
or, in the (r) representation, 



M a Mp = {Tp - r a )F'(0) + - (r| - r Q 2 ) F"(0) + ... 

and 

F'(r^ = F 7 (0) + r a F T iO) + ... 

It follows 

F 7 ^ A Mjtp = \(r a - Tpf F^JO) A F^(0) + ... 
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When r a and Tp tends to 0, the plane (M a , M a Mp, F'(tJ)) tends to the oscu 



i 2 F'(r a ) A M a Mg , , ,. . 

latory plane to (1 ) at G because the vector -r — - admits the limit 

{T a - Tp) 2 



F'(0)AF"(0). 

Theorem. The parametric representation of curve T in Frenet frame at C is 
x = t , y = t 2 , z = t 4 ; rG t>(0) . 

Proof. 

Let C an ordinary point of (T). In a neighborhood of C, we have 

x = r,?/ = r 2 ,z = r 3 ; r€ u(0) (24) 

such that t = corresponds to C. Let M Q and Mg two points of T with 
r a Tp < 0. Eqs. (fJU) lead to 



F'{r a ) A M a Mf : 



and by changing rg with r Q 



F'fo) A MpM 



( „3 _|_ o_ _2 _2_ \ 

2"^ - r| - r a Tp 

Tp - T a 



/_ t 3 + 2t/3T 2_ TqT 2\ 

2t| - r 2 - r^r a 



(25) 



V 



The two vectors in Eqs. 



T a - Tp 
are parallel if and only if 



(26) 



rl + t} - 2r a rp = 



'a ~ '/3 



Ta T n T, 



■ a Tp - TpT a 



0. 



(27) 



Equation fpZTjb yields r Q = rg and then (J27j) 2 is verified. This means that any 
relation Tp = <f)(T a ) different from identity is not allowed and r a Tp < is not 
possible. Then, the critical point C cannot be an ordinary point of (r). Let 
us now consider the simplest singular case for T at C 



2 4 
X = T , y = T , Z — T ■ 



The parallelism conditions for vectors F'(r a ) A M a Mp and F'(rg) A MpM a are 
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equivalent to 

-2r Q 4 - 2r 3 a r p + 2t 2 t 2 + 2r a rj = 2r| + 2rjr a - 2r|rJ - 2r p r 3 a 
3^ - r 2 a rp - r a r 2 - r| = -3rJ + r 2 r a + r^r 2 + . 

Equation ([2"51) 9 implies 2(r Q + r / g)(r Q , — r^) 2 = which has, for r a ^ rg, the 

unique possible solution r a = — tq and then (pZBjb is satisfied. 

The equations of the projection (T f ) of the curve (T) on the plane y = r 2 are 

4 

x = r, z = r 

and the straight line M a Mp remains bitangent to (T'). 

Consequently, in the plane y = r 2 , the main part of the equation of (T') is 
written in the form 

z = H{r a ) {x - r a f (x + r a f + tI with H(0) ^ 0. 

In adimensional form, H(0) = 1 and the main part of the equation of surface 
E generated by this curve in a neighborhood of the critical point C, when 
point C is taken as origin of the axes and Cxy is the osculatory plane to S at 
C, is 

z = (x 2 - y) 2 + y 2 . 

The dimensions of x and y are different. In an adimensional system of axes, 
the equation writes 

z = A \(B 2 x 2 -y) 2 + y 2 
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